Spectrum of Light in a Quantum Fluctuating Periodic Structure 
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We address the general problem of the excitation spectrum for light coupled to scatterers having 
quantum fluctuating positions around the sites of a periodic lattice. In addition to providing an 
imaginary part to the spectrum, we show that these quantum fluctuations affect the real part of the 
spectrum, in a way that we determine analytically. Our predictions may be observed with ultracold 
atoms in an optical lattice, onaJ = O^J' = l narrow atomic transition. As a side result, we 
resolve a controversy for the occurrence of a spectral gap in a fee lattice. 

PACS numbers: 42.50.Ct, 67.85.-d, 71.36. +c 



The investigation of light propagation in a periodic 
structure is a fundamental problem in condensed mat- 
ter physics P, 0] , ranging from the physics of photonic 
crystals Q to X ray Q or even 7 ray 0, 0] scatter- 
ing by a crystal. It recently gained a renewed interest 
0, i,S [13, El, El thanks to the possibility of produc- 
ing artificial periodic structures of quantum dots [9| or of 
atoms [isj , which have a much larger spatial period than 
natural crystals, allowing a measurement with a laser of 
the excitation spectrum over the whole Brillouin zone. 

In reality, strictly periodic structures do not exist: the 
quantum (if not thermal) fluctuations of the positions 
of the scatterers in the structure are unavoidable. From 
general grounds, we know that these fluctuations affect 
the spectrum of light in two ways. First, they introduce 
a dissipative component: the elementary excitation spec- 
trum acquires an imaginary part, a well known effect of 
phonon coupling in condensed matter physics. Second, 
they introduce a reactive component, modifying the real 
part of the excitation spectrum. Whereas the fluctua- 
tions of the scatterers positions have indeed been taken 
into account in simplified models [ll, the only explicit 
prediction for the corresponding modification to the ex- 
citation spectrum was given in [5| . 

The problem however is not closed yet: in the limit 
of vanishing position fluctuations, Eq.(4.9) of [1] for the 
real part of the spectrum does not reduce to the predic- 
tion given in for fixed scatterers, a prediction that also 
does not coincide with the one of [Tj ■ This problem is not 
only formal, it may soon be addressed in current exper- 
iments with narrow line cold atoms trapped in optical 
lattices 1^, where measurements of the spectrum with 
a good precision may be performed. Here we provide a 
conclusive analytical answer to this problem. 
Model: Although our method to come is general, we as- 
sume for concreteness that the scatterers are atoms cou- 
pled to the electromagnetic field on an electronic tran- 
sition between the ground state g of spin and an ex- 
cited state e of spin 1. The atoms are trapped at the 
nodes of an optical lattice, with nowhere more than one 
In the deep-lattice limit, tunneling 
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atom per site 

is negligible and the i*"^ atom is assumed to be harmon- 



ically trapped around lattice site R^, with the potential 



^ho 



(fi — R;) /2. Here m is the atomic mass, 



is the position operator of atom i, and Who is the atomic 
oscillation frequency. In this regime, the fluctuations of 
the atomic positions are purely on-site and uncorrelated, 
contrarily to the case of phonons in a crystal The 
Hamiltonian of our system [15[ may be split in the non- 
interacting term Hq and the atom-field dipolar coupling 
V, H ^ Ho + V, with 



Ho 



N 

E 



g- + c/,(f,) + ^/i^oN 



/ d'^ k hck a|j.g flke . 



(1) 



Here N is the number of atoms, loq is the bare atomic res- 
onance frequency, the sum over the three directions 
of space X, y, z accounts for the three-fold degeneracy of 
e, I? is the three-dimensional Fourier space truncated by 
a cut-off fc < fc^, and the annihilation and creation oper- 
ators obey usual bosonic commutation relations such as 
= (5£.£'(5(k — k'), where e and k are the photon 



polarization and wavevector. The coupling V is 



N 



V : 



E±(f, 



(2) 



where Di,a = d\i : ea){i : g\ +h.c. is the component along 
direction a of the dipole operator of the i"^ atom, 
proportional to the atomic dipole moment d, Ej_(r) = 
/■p ci'^fc X]e±k [^k£ akE e*'^ '" -f- h.c] is the transverse elec- 
tric field operator, and £k = i(27r)"^/^[fi.A:c/(2eo)]^^^- 
Method: We treat the coupling V to second order of per- 
turbation theory, to calculate atom-field elementary exci- 
tations mainly of atomic nature. In practice for a lattice 
spacing ~ 1 /fco, where fco = luq /c is the resonant wavevec- 
tor, this requires F <C ojho, where F = d'^ / {STreofi) is 
the free space atomic spontaneous emission rate. The 
energies of the system up to second order in V are eigen- 
values of the effective Hamiltonian 

Q 



H^s = PHP + PVQ 



Eio)Q-QHoQ 



QVP, (3) 
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where E'^"^ is an eigenenergy of Hq, P projects orthog- 
onally onto the corresponding eigenspace of Hq, and 
Q = I — P. We now obtain the excitation energies of 
the system as the difference between excited states ener- 
gies and the ground state energy. 

For the perturbative calculation of the ground state 
energy Eg, we take for P the projector over the state 
with all the atoms in the electronic and motional ground 
states, and the field in vacuum, so that E, 
We obtain Ea 



(0) 



Ncg, with 



-huj] 



ho 



rf2 



(fik 



hck 



V 



(2tt)^ hujQ + hck 



(4) 



For the excited state energies, we take for P the projector 
Pf, over all the states ||i : Bq), where ||i : Ca) represents 
the atom i in the electronic state 6^, the iV — 1 other 
atoms in the ground electronic state, all the atoms in the 
motional ground state, and the field in vacuum. Then 
P projects over a subspace of dimension 3iV. Since this 
excited subspace is coupled by V to the continuous part 
of the spectrum of _ffo, one has to replace i?*-"-* by E^^ + 



iO+ in (l3|), with Ee~Eg+ hojo, which gives rise to a 
complex excited state energy Eg, eigenvalue of 

i?,^ff=[(iV-l)6,+ee]Pe 

+ ■ ea)(j : ef3\\. (5) 



Here eg is the complex energy of a single atom [17 1 

3^ £ f (Pk hck 

2 3eo Jv {^T^y "^0 + «0+ - hck 



(6) 



To obtain the excitation energies of the system we sub- 



tract the ground state energy Eg from 



In this sub- 



traction, the dangerous terms proportional to the num- 
ber of atoms N disappear and the interatomic distance 
independent term gives the excitation energy for a single 
atom, that we split in a real part and an imaginary part: 



huJA 



hv 



(7) 



As expected, F is the free space spontaneous emission 
rate, and the effective atomic resonance frequency loa 
deviates from wo by Lamb shift type terms. 

Most interesting are the position dependent terms 
in ([5]), which contain the effective coupling amplitude 
3a/3(r) for the transfer of the atomic excitation in be- 
tween two different sites separated by r. This cou- 
pling amplitude appears as the inverse Fourier transform 
5Q/3(r) = /p[d^fc/(27r)^]e'*' '"gQ,;3(k) of the function 



5a/3(k) = 



37r?iF k'^Sap — kakp 



-k~at 



h.3, 



k^ — k'^ 



iO+ 



(8) 



The effect of the quantum fluctuations of the atomic 
positions on the intersitc coupling here enters through 



the size of the harmonic oscillator ground state Oho = 

[ft/(2mc.ho)]^/'. 

The integral defining gapiv) is cut at large k by the 
Gaussian factor of momentum width 1 /oho k^ ■ Hence 
we can evaluate the coupling amplitude Qapi^i — R-j) by 
extending the integral over k to the whole space. In 
this limit the coupling amplitude is the average, over the 
harmonic oscillator ground state probability distributions 
of Vi and r^, of the function gap^TCi — fj), where [lit 



fa/3(r) 



■ikor 



(9) 



is proportional to the component along a of the classical 
electric field radiated by a point-like dipolc of frequency 
luq oriented along direction /3 [lit ■ ^ crucial consequence 
of the average is that, whereas gap{j^) has the electro- 
static divergence at the origin, the function g^piv) 
is regular even in r = 0, with a value 



_ , , hT ^ 
gai3(0) = —Oap 



Erfi (fcoflho) ~i 1 + 2(A;oaho)^ 



27ri/2(fcoaho)3 



(10) 

where Erfi (x) = 27r ^^'^ dy exp(?/^) is the imaginary 
error function. For large r, the average over the atomic 
motion docs not suppress the long range nature of the 
radiated dipolar field cx exp(iA:o'')/''- From (jl7p one has 
indeed gapiic) — 5a;3(r) exp(— /CqO^q) remarkably as soon 
as r » Oho- 

Periodic case: We now take the limit N +00 with one 
atom per lattice site, realizing for the light field a periodic 
potential, here an arbitrary Bravais lattice. To diagonal- 
ize H^g we rely on Bloch theorem: The eigenvectors ji/^q) 
depend on the lattice site position as 



(11) 



where the Bloch vector q is chosen in the first Brillouin 
zone of the lattice. I.e. the "dipole" carried by atom i 
differs from the one d carried by the atom in R = 
by a global phase factor. Thus the infinite dimension 
eigenvalue problem on the excitation spectrum £q. 



[H:g~EgP,]\^JJ^)^e^\lP^), 



(12) 



reduces to the diagonalization of the 3x3 matrix M, 



Md = Eqd, with matrix elements 
Map = (hujA ~ ) Sap + ^ 



apiR)e'''-^, (13) 



where the sum runs over the lattice L excluding the ori- 
gin. By adding and subtracting gap{0) and using Pois- 
son's summation formula we convert this sum into a sum 
over the reciprocal lattice RL: 



Map = huJA 



■^r\ r 1 
I— I Oap-gap[0) + 



Vl ^ 



5a/3(K-q) 
(14) 
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where Vl is the primitive unit cell volume of the lattice. 
Imaginary part: Since each term of the sum over K in 
is real, the imaginary part of the excitation spectrum 
can be calculated explicitly: 

hT 



Im Ec, 



1 



(15) 



As already mentioned, this non-zero value is a direct con- 
sequence of the fluctuating scatterers positions. This is 
why an expression analogous, but not equal to (jisp was 
derived in Q in the different context of 7 ray nuclear 
scattering in a crystal. 

The non-zero imaginary part of Sq is indeed due to 
the decay of the system out of the subspace where all 
the atoms are in their motional ground state. When 
an excited atom Bq in its motional ground state |0)ho 
spontaneously emits a photon of polarization e and mo- 
mentum fcon, |n| = 1, its probability density to fall in 
g with an excited motional state is [18| (3r/2)|eQp[l — 
|ho(0| exp(— ifcon • f)|0)hop], which after sum over e _L n 
and average over the direction n, exactly gives the de- 
cay rate — 2Imeq//i. This process conserves the quasi- 
momentum q. If the emitted photon carries away the 
quasi- momentum Qph, the resulting g atom in the mo- 
tional excited state is coherently delocalized over the 
whole lattice, with a probability amplitude oc e'l^^f^ of 
being in site R, and a quasi-momcntum Qat = q — qph- 
Since qph belongs to a continuum, this spontaneous emis- 
sion process opens up a continuum of final states, hence 
the possibility to have for a fixed q a continuous spec- 
trum for H and a complex £q. Experimentally, to obtain 
long lived elementary excitations, one may operate in the 
so-called Lamb-Dicke regime, koUho ^ 1, where the loss 
rate — 2Imeq/fi, ~ r(/cottho)^ is much smaller than T. 
Real part: In we replace ga/siO) and gap by their 



explicit expressions (|10p and ([8j) . Then Re Sq — hujA is an 
eigenvalue of the 3x3 real symmetric matrix 



SMaf: 



hv 



Saf: 



1 + 2(A:oaho)^ 



37r;ir 



27ri/2(fcoaho)3 



Erfi (fcoflho) 



E 



X'2 



(16) 



where K' = K — q [19|. Eq. lfTHl) is useful for a numeri- 
cal calculation of the spectrum, see Fig[TJ Remarkably, 
one can even derive analytically the dependence of the 
spectrum with Oho from (|13p : one multiplies gafsiR) by 
e'^o^ho and one takes the derivative with respect to a^^. 
In the Fourier representation ([8|), this pulls out a factor 
/cq — fc^ which exactly cancels the denominator. The re- 
sulting Fourier integral is now essentially Gaussian and 
can be directly calculated: 



i/3(r) = 



Jo' 



8^1/2 (fc„ah< 



-Ja/jAr -I- dr^dris) 6 



(17) 
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FIG. 1: Spectrum of light in a quantum fluctuating periodic 
atomic structure, (a) For a simple cubic lattice, with a lattice 
constant a (Vl =1"^), comparison of the analytical prediction 
(|19p (solid lines) with the numerical solution of (|16p (sym- 
bols), for two values of koa and of the Bloch vector (x : point 
R and o: point F of the first Brillouin zone), (b) and (c): 
Real part of the spectrum as a function of the Bloch vector 
along the standard irreducible path in the first Brillouin zone, 
for (b) a simple cubic lattice, and for (c) a face-centered cubic 
lattice with a lattice constant 2a (Vl = 2a"'); here koa — 2 
and fcoctho ~ l/\/30, which leads to aho/i — 0.09. In both 
cases the band structure is gapless. 



From (|13p one then obtains a sum excluding the origin 

7 2 



SMa 



E 

RGL* 



Mq/3(R) 



(18) 



It is apparent that each term of the sum is exponentially 
small since, in the deep lattice limit, the harmonic length 
flho is much smaller than the lattice spacing. Keeping 
these terms is actually beyond accuracy of our Hamil- 
tonian (see [iBl)- One can thus set the right-hand side 
of (dll) to zero, which leads to SMap = e-'^o^L SM^^^, 
where SM^p is independent of a^o and is simply the limit 
of SMap when Oho 0. This shows that the real part of 
the excitation spectrum is a Gaussian function of aho^ 



RcEq — ftLUA 



e 



hujA) 



(19) 



where ; 



is the limit of the excitation energy for Oho — > 0. 
Eg. pop is the main result of this work, it gives in a very 
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explicit way the influence of the fluctuations of the atomic 
positions on the real part of the light spectrum. It is in 
excellent agreement with the numerics, see FiglT^. 
Band structures: We illustrate by numerically cal- 
culating the band structure in the two lattice geometries 
of 0, M- III FigH we show Re £q — fkuA as a function of 
the Bloch vector along the standard irreducible path in 
the first Brillouin zone, for the simple cubic (Fig[T}j) and 
for the face-centered cubic (fee) (Fig[Tt) lattices. The 
figure reveals the lack of an omnidirectional gap, which 
confirms the prediction of Q against the one of Q for 
fixed atomic positions [2^. We also explored the bee 
and several generic less symmetric Bravais lattices with- 
out finding an omnidirectional gap. On the contrary, 
the non-Bravais atomic diamond lattice may support an 
omnidirectional gap Note that, according to 

changing Oho amounts to a mere rescaling of the vertical 
axis of FiglT]D,c and cannot open or close an energy gap. 
The situation may be different for anisotropic microtraps. 
Experimental issues: The perturbative regime F <^ Who 
considered here was not usual in atomic lattice experi- 
ments. It is now available in experiments using very nar- 
row transitions for atomic clock purposes: e.g. ^^Sr was 
recently trapped in a deep 3D lattice, with Oho ^ 0.05a ~ 
20nm [IJ], and it has a narrow line 5s^ — > 5s5p ^Pi 
realizing the needed J = ^ J' = 1 transition, with 
F ^ O.OSwho- To produce and spectroscopically probe 
elementary excitations with q ^ kg, one cannot use the 
direct g — > e coupling with resonant light, but one can 
use an indirect Raman coupling [2ll |. 

We also assumed that g and the three sublevels of e 
experience the same trapping potential. This is the case 
in Wigner ion crystals [22| , where e and g experience the 
same Coulomb shift. For neutral atoms, the lattice po- 
tential is a lightshift, that may deviate from one of the 
two assumptions of ([1]), (i) e experiences a scalar light- 
shift, and (ii) the lightshifts of e and g arc equal. For 
a lattice obtained by incoherent superposition of laser 
standing waves along x, y, z linearly polarized along y, 
z and X respectively, violation of (i) breaks the harmonic 
oscillator isotropy: the sublevel (resp. ey,ez) has an 
oscillator length rjeCiho along z (resp. x, y) different from 
the one aho along the other two directions. This does not 
break the three-fold degeneracy of the motional ground 
state in e but it reduces the overlap between the mo- 
tional ground state of Cq and that of g. Hence after 
spontaneous emission, even if one neglects the atom re- 
coil (/cQfiho —>- 0), the atom in g can populate an excited 
motional state, giving a non-zero decay rate to the ele- 
mentary excitations. Similarly, violation of condition (ii) 
leads to an oscillator length in g equal to riga^o, Vg 
which also increases the decay rate. For fcoiho 0, com- 
bining both violations gives 



2 V ^ [VeVg + iVeJlg)^^]{Vg + Vg ) ^ 

E.g. if one has achieved r^g ~ 1 at the expense of having 



an optical lattice depth in g twice as small/large as in e, 
one still finds a small decay rate ~ 0.04F. 

Conclusion: Quantum fluctuations of the positions of the 
scatterers in a periodic structure very generally give rise 
to an imaginary part in the spectrum of light and affect 
its real part. For scatterers tightly trapped in a periodic 
potential, we derived an expression for this spectrum. We 
showed that, amazingly, its dependence on the amplitude 
of the fluctuations of thepositions is a Gaussian, not only 
for the imaginary part [5|] but also for the real part. This 
effect on the real part can be large and may be observed 
in recent atomic lattice clock experiments. An intriguing 
perspective is the extension of this work to the disordered 
case and to localized states of light. 
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